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A mapping of avalanches occurring in the zero-temperature random-field Ising model (zt-RFIM) 
to life-periods of a population experiencing immigration is established. Such a mapping allows the 
microscopic criteria for occurrence of an infinite avalanche in a g-regular graph to be determined. 
Based on these criteria, we explain why an infinite avalanche can occur only for g > 3. A key factor 
for an avalanche to become infinite is that it interacts in an optimal way with previously flipped 
spins. The generating function techniques developed for branching processes are applied to obtain 
expressions analytically for the duration, pulse-shapes and power spectrum of the avalanches. The 
results show that only very long avalanches exhibit a significant degree of universality. 

PACS numbers: 75.10.Nr, 64.60.aq, 75.60.Jk 



Systems exhibiting avalanche response when driven ex- 
ternally are ubiquitous and have been widely studied in 
recent decades [l[ . Examples of such phenomena include 
magnetisation reversal 0, fracture and earthquakes [H or 
collective opinion shifts [J]. Broadly speaking, this kind 
of systems can either pop (all the avalanches are of sim- 
ilar size), snap (there is at least one characteristic large 
event) or crackle (avalanches have a broad range of sizes). 
The competition between spatial interactions promoting 
snap behaviour and heterogeneity that favour popping 
has been proposed as a key factor dictating whether a 
system will snap, pop or crackle [ij. However, the mech- 
anisms which determine the nature of such avalanches 
are still poorly understood and this affects the ability to 
make predictions. For instance, it is often a challenging 
task to predict the occurrence of large avalanches (e.g. 
a devastating earthquake or an abrupt collective change 
of opinion) by observing the evolution of systems over a 
short period of time. 

In this Letter, we investigate the link between micro- 
scopic dynamics and the type of global avalanche re- 
sponse within the framework of the zero-temperature 
random- field Ising model (zt-RFIM) P, This proto- 
type model, originally proposed for disordered magnetic 
materials, provides a unified description of the three com- 
mon avalanche types in various systems. We demonstrate 
that the time-dependent behaviour of avalanches in the 
zt-RFIM defined on locally tree-like topologies (e.g. a 
q-regular graph, where q is the coordination number) is 
described by a branching process (BP) which is a well- 
established model for population dynamics [1] . The map- 
ping provides the criteria for the appearance of an infi- 
nite avalanche in the system. In particular, we give an 
explanation for the long-standing question 0, Q of why 
infinite avalanches can occur in a Bethe lattice or a q- 
regular graph with g > 3, but do not occur for q < 3. 
In addition, the proposed framework leads to exact so- 
lutions for experimentally important quantities such as 



the distribution of avalanche durations, their mean pulse 
shape and mean power spectrum. 

The Model - The zt-RFIM consists of a set of N spin 
variables = ±1 placed on the nodes of a network with 
links representing interaction between pairs of spins. The 
system is described by the following Hamiltonian, 
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where the first summation is taken over all pairs of inter- 
acting sites, (j,j) and interaction strength J ~ 1. Het- 
erogeneity is introduced by quenched and independent 
local random fields, hi, distributed with the probabil- 
ity density function (p.d.f.) p{hi), which is continuous 
and non-zero everywhere and characterised by a zero 
mean and typical variance A^. The system is driven 
by the external field, H{t), which is assumed to in- 
crease monotonically and adiabatically from — oo (when 
all spins are in the down-state, Si = —1) to -l-oo. Spins 
flip to the up-state (si = -1-1) at a rate F according to 
single spin-flip dynamics [^, i.e. they align with their 
local fields, fi{t) = H{t) + hi + J{2ni{t) — q), where 
ni{t) = J2j/i i^j + 1) /2 is the number of spins j which 
neighbour spin i and are in the up-state. With this dy- 
namics, the system exhibits intermittent response with 
popping being observed for any network topology at large 
disorder. An infinite avalanche occurs at low disorder for 
certain network topologies: a complete graph (mean-field 
model 3|), a Bethe lattice or a g- regular graph with coor- 
dination number q > 3_|9(, and on hypcrcubic lattices in 



dimensions d > 3 



For these topologies, pop and 



snap regimes are separated by a disorder-induced contin- 
uous phase transition located at a topology-dependent 
critical degree of disorder, A = Ac, where crackling oc- 
curs. In contrast, no disorder-induced transition is ob- 
served for g-regular graphs with g < 3 which exhibit 
popping for any non-zero disorder 
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Transformation dynamics of the zt-RFIM. - In order to 
establish the time-dependent behaviour of the zt-RFIM 
and map it to a BP, it is convenient to introduce three 
states for spins: stable (S), unstable (U) and flipped (F). 
As the external field is driven from —oo to oo, each spin 
passes consecutively through these states, S— >U^'F (see 
an example in Fig. [1]). When spin i is in the initial state, 
S, the field is negative (/i < 0) and the spin is down 
{si = — 1). At some time, the local field at spin i becomes 
positive {fi > 0) and the spin goes to state U meaning it 
becomes unstable but is still in the down-state. Finally, 
at a later time, spin i flips and moves into the flipped 
state F characterised by > and Si = +1. The first 
transition, S— >U, is mediated by the local field. As such 
it can be caused by one of two mechanisms: (i) field- 
induced, i.e. due to the increase in external field H{t) 
or (ii) flip-induced, caused by the flip of a neighbouring 
spin and associated increase in ni{t). The second tran- 
sition, U— >-F, occurs as a Poisson process at rate F. As 
a consequence of the second transition, some number 
of the neighbours of spin i move from state S->U. This 
dynamical behaviour can be viewed as a BP describing 
an evolving population in which two types of event can 
take place: immigration of new individuals and reproduc- 
tion followed by immediate death Q. Indeed, the spins 
in state U form a population of individuals that are cre- 
ated spontaneously (immigration) by the fleld-induced 
mechanism and can produce new U spins (reproduction) 
when passing to the state F (death) according to the flip- 
induce mechanism. Both immigration and reproduction 
are stochastic, and follow Poisson dynamical rules, with 
immigration occurring at a rate and reproduction at 
a rate F for each individual. In the example shown in 
Fig. [Jl an immigration event occurs at time t2 followed 
by the reproduction events U U 2U. In a BP, 
the number of offspring produced by an individual is a 
random number ^ > 0. If the expectation value E[^] < 1, 
then the population size will be non-zero only for cer- 
tain "life-periods" of short duration At the end of 
each life-period the population becomes extinct, with the 
next life-period being started by an immigration event. 
Such life-periods correspond to popping in the zt-RFIM 
If E[^] > 1, the population size can grow exponentially, 
analogous to an infinite avalanche of spin flips in the zt- 
RFIM, causing the system to snap. 

Infinite avalanche in a q-regular graph.- From popu- 
lation dynamics it is clear that an infinite avalanche can 
occur only when the mean number of offspring E[^] > 1. 
However, applying this criterion to the zt-RFIM is chal- 
lenging because E[^] is a function of time which depends 
on H{t) and is restricted by the time-dependent avail- 
ability of S spins in the neighbourhood of U spins. This 
contrasts with BPs proposed for description of avalanches 
in critical stationary states with constant E[^] [l^. [Tsj. 
For the zt-RFIM, E[^] = at the initial moment of time. 



t = 



t = t-i 



U 



t = t2 



U 



U 



u 



FIG. 1. Example of an avalanche on part of a 4-regular graph. 
At time ti, the configuration consists of only S (stable) and F 
(flipped) spins. As the external field increases, a spin becomes 
unstable (U) at time t2 > t\. This represents an immigration 
event for the population of U sites. At time ts > t2, the U spin 
flips causing the central site, i, to undergo a 5* — >■ [/ transition. 
The central site, i, had one flipped neighbour (along link 4) 
prior to time ts, i.e. ni{t < t^) = 1, and it became unstable, 
so that C^i is equal to the number of neighbouring sites in 
state S, i.e. (i = 2. In contrast, C^i would have been zero if 
the central site had remained stable. At time t4, the central 
site flips and this leads to 2 new U offspring. 



when H ~ —oo, and an infinite avalanche only occurs if 
it increases to become greater than one. The deriva- 
tive of E[^] with respect to time consists of two con- 
tributions associated, respectively, with flip- and field- 
induced changes: E[^] = Ait)Cov[Ci,n,]+Bit) Here, 
Cov[G, Hi] = E^[n^{Q ~ EP[Cj])]. E^[. . .] denotes a mean 
over all stable neighbours i of flipping spins, while ni(t) 
is the number of neighbours of each spin i which have 
already flipped (pre-flipped spins) and Q is equal to the 
number of stable neighbours, q — I — ni{t), of spin i if 
it also flips and 0, otherwise (see Fig. [T]). It can be 
shown that E^[Ci] = E[^] The function B{t) > 

remains finite at all times. In contrast, A{t) > diverges 
as A{t) oc (I — E[^])~^, meaning that the fiip-induced 
mechanism is the main contribution to E[^] when E[^] 
approaches 1 from below. In this case, A{t) ^ 1 and the 
system can only become supercritical if E[^] > 0, mean- 
ing that the following condition must necessarily hold: 



lim CovfCi, nA > 



(2) 



This criterion constitutes one of the main analytical find- 
ings of this Letter. 

For a g = 3 lattice, Cov[C^,n,] = {I -E^[Q])P^ {1, 1) - 
E^[C,]{P^{0,l) + 2P^{0,2)), where P^{k^) is the prob- 
ability that Q = k and n.i = m (see [16|). The only 
positive contribution (the first one), tends to zero as 
E^[Ci] 1", while the negative contributions do not, 
meaning that Cov[Ci,ri,i] becomes negative and condi- 
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tion ([2]) is not satisfied. The reason for this is that the 
interaction of an avalanche with pre-flipped spins (occur- 
ring at spins where Ui > 1) must cause the avalanche to 
branch (i.e. Q must be greater than one for this spin; see 
the evolution of the state of the central spin in Fig. [T]), 
thus making the BP super-critical. This type of branch- 
ing requires at least 3 open paths (see paths 1, 2 and 3 in 
Fig- HI) from a flipping spin (one incoming (path 1) and at 
least two outgoing (paths 2 and 3)) in addition to at least 
one path closed by the presence of the pre-flipped spins 
(path 4). This means that an infinite avalanche driven 
by interaction with pre-flipped spins is only possible for 
q> 4= [1 (closed path) +1 (incoming) +2 (outgoing)]. 

Numerical support for condition ([2]) is presented in 
Fig. [21 The avalanches below (a,b,c) and above (d,e,f) 
criticality are analysed for a 4-regular graph. The infi- 
nite avalanche occurs stochastically after E^[Ci] (the rate 
of branching) passes 1 (see (a)), at which time there is 
a non-zero probability of any given avalanche being in- 
finite. The mean positive derivative of E^[^j;] prior to 
the infinite avalanche is due to the fact that Cov[C,;, n^] is 
positive when E^[Ci] 1~ (see (a) and (b)). The infinite 
avalanche causes the covariance to reduce and eventually 
become negative, which, in turn leads to E[^] dropping 
below 1 and the avalanche terminating (see (c)). Above 
the critical degree of disorder (A > Ac), the branch- 
ing rate grows initially and, despite approaching 1 it is 
always < 1 (see (d)). This is consistent with the fact 
that condition ([2]) is not satisfied since Cov[Ci,ni] be- 
comes negative before [Q] reaches its maximum value 
(cf. (d-f)). These results suggest that monitoring the 
evolution of E^[Ci] and Cov[(i,ni] over relatively short 
time intervals prior to observing a large avalanche could 
be used to predict its occurrence. 

Time- dependent properties of avalanches.- Using the 
known generating function techniques for BPs Q, it is 
possible to derive analytical expressions for the power 
spectrum, the distribution of avalanche durations and 
pulse shapes. These experimentally relevant quanti- 
ties 17h20|| can be defined in terms of the rate of change 
of magnetisation, d(m)/di (where (m) = {J2i^i))^ 
which is proportional to the experimentally measured 
voltage d{m)/dt cx Vit) = VoJ^i^it - ti) (with Vq 
being a coefficient of proportionality) where each S- 
function corresponds to the flip of a spin at time ti. 
The power spectrum of avalanches, defined as the mean 
(over all avalanches occurring near external field H) 
square of the Fourier transform of the voltage {S{lu)) = 

can be found exactly for a 
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FIG. 2. The rate of branching E^[0] (a,d) and Cov[Ci,ni] 
(b,e) vs time t. Panels (c) and (f) show Cov[(^i,ni] vs 
for data presented in (a)-(b) and (d)-(e), respectively. Data 
corresponds to a 4-regular graph of A'^ = 10® spins, with nor- 
mally distributed disorder below A = 1.65 (a,b,c), and above 
A = 2.5 (d,e,f), critical Ac ~ 1.78215. The external field was 
swept from — oo to -l-cxa at a rate AH /At = 2 x 10"® with 
H{t = Q) = Q and stochastic flips occurred at a rate F = 1. 



functional form as that found within the mean-field ap- 
proach 17, 1^, with the additional constant factor Sq. 
At critical disorder A = Ac and near the critical ex- 
ternal field He, the value of Si oc \H — Hc\ ^^"^ and 
the frequency-independent term in the denominator ap- 



H,f\ i.e. {S{uj)) cx 
This divergence in 



where the frequency-independent factors, Sq and Si are 
given in Ref. [lij . The power spectrum has the same 



proaches zero according to \H 
Si{H)uj-^ (see Fig. Hb)) ^ ' 
S{uj) can be integrated out finding the average of S{uj) for 
all avalanches occurring in a single magnetisation rever- 
sal, i.e. (5'int(a;)) = /^^_^(5(a;))i?Mjy. This quantity 
also follows a power law at criticality, (S'int(w)) cx w"^, 
with = 3/2 (see inset of Fig. [HJa)). This exponent can 
be compared with the experimentally observed exponent 
in the range 6* ~ 1.6 - 1.7 [H,!!!. 

Let the first spin fiip in a single avalanche occur at 
t = and the last one at t = T. The duration T of 
the avalanche obeys a field-dependent p.d.f. pT given by 
(see [l^ for detail), 

PT =a,a(/T(0))9T/T(0) , (4) 

where a{x) = X]|i=o *'^^'^^(^^) ^® generating function 
for the number, of spins destabilised by the first spin 
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flip in the avalanche. The generating function frix) = 
E[a;'^^] describes the number, Xt^ of spins which are 
unstable &it ~T due to the influence of a single unstable 
spin at t = 0. Away from criticality, pT oc exp(— T(l — 
E[^]), and at the critical point, pr oc (se e FigJHTb)), 
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10 10 10 



' (see 
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which is the same as found in mean-field 

The pulse shape {V{t))T is the mean voltage measured 
for all avalanches of duration T happening near a partic- 
ular value of external field, and is given by (see [16| for 
detail) , 

{y{t))T^v,T ('%#^a./.(/T-.(o)) 



d.Mh-tm 

+VoSit) + VoSiT - t) 



(5) 



where f{x) = X]|=o •^^-^(^) generating function 

for the number, ^, of spins destabilised by the flip of 
any spin in the avalanche except the first one. At crit- 
icality, the shape is similar to that for a simple random 
walk (sec Fig.[3fc)). while away from criticality, the pulse 
shape flattens for large T, similar to the effect found 
for a ID random walk in a parabolic potential well [26 1 
(see Fig. |3jd)). At A = Ac, the pulse shape is signif- 
icantly asymmetric for T < 100 and a scaling hypoth- 
esis 0, [S H li^l of the form {V{t))T = TW{t/T) is 
only appropriate for durations T > 100 The distri- 
bution pt also deviates significantly from the power-law 
behaviour for short avalanches (Fig. [UJb)). Therefore, 
experimental studies might not reliably detect power-law 
scaling of dynamical properties unless they focus on very 
long avalanches, with the drawback that they are much 
less frequent. 

Conclusions.- The mechanisms of evolution of 
avalanches in the zero-temperature random field Ising 
model on a q-regular graph have been explained in terms 
of a mapping to a branching process describing an evolv- 
ing population with immigration. The global avalanche 
response has been linked to local quantities such as the 
mean number E[^] of spins becoming unstable after the 
flip of a neighbour and the interaction between the prop- 
agating front of an avalanche and spins flipped in previ- 
ous avalanches (formally accounted by Cov[Q, rii]). From 
this, we have answered the question as to why the jump 
in magnetisation does not occur in Bethe lattices with 
g < 3 p, S| • The presented analysis suggests that moni- 
toring the evolution of appropriate local quantities can be 
used to predict the occurrence of infinite avalanches. Our 
results are based on a relatively simple framework but 
they serve as a solid mathematical foundation for future 
work aiming at elucidating the evolutionary mechanisms 
of avalanches in more realistic systems. 

TPH acknowledges the financial support of the EP- 
SRC. 
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FIG. 3. Panels (a) and (b) show the integrated power spec- 
tra, (Sint(<^)) and the distribution of avalanche durations, pr, 
respectively, for a 4-regular graph with normally distributed 
disorder, A = Ac (critical disorder; solid lines and circles), 
A = 2.0 (dashed lines and squares), A — 2.5 (dot-dashed 
lines and diamonds) and A = 3.0 (double dot-dashed lines 
and triangles). The inset of panel (a) shows the power spectra, 
{S{ij)} evaluated for the same parameters as for (a) and (b). 
In both (b) and the inset of (a), H — 1.0. In (a) and its inset, 
the dotted lines indicate S'int(i^) oc uj~^^^ and S{ui) oc o;"^, 
respectively, corresponding to the asymptotic behaviour for 
uj <^ T near critical disorder. The solid line in (b) tends to 
the limit pr oc for T ^ P"^, shown by the dotted line. 
Panels (c) and (d) show the shape of avalanches {V{t))T vs 
t/T for A = Ac and A = 2.5, respectively, H = 1.0 and 
T = 10 (solid lines and circles), T = 20 (dashed lines and 
squares) and T = 30 (dot-dashed lines and diamonds). Ana- 
lytical calculations are represented by lines. Numerical data 
(represented by symbols) for durations and spectra were av- 
eraged over avalanches occurring when, 0.998 < H < 1.002 
(inset of (a) and (b)) and 0.98 < H < 1.02 ((c) and (d)) for 
10^ realisations of 10^ spins and P = Vb = 1. 
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I. MAPPING OF ZT-RFIM TO A BRANCHING PROCESS 

The mapping between the zt-RFIM and the branching process (BP) described in the main text is 
quantified by the following values: (i) the rate, R^, of spontaneous transitions S— )-U in field-induced 
changes, i.e. the rate of immigration in the BP; (ii) the probability, P^^ , that a stable spin becomes 
unstable, S— t-U, when its neighbour changes from U to F (a flip-induced change describing reproduction); 
(iii) the distributions, P^{Ci) and P^{Ci), of the numbers of neighbours, Q, of site i which are in state S 
immediately after site i moves from S— 7>U according to either a field- or fiip-induced change, respectively 
and (iv) the distributions P^{^^) and P^(^^) of the numbers of offspring of sites produced by field- and 
fiip-induced changes, respectively. In what follows and in the main text, the superscript R in is 
suppressed, and P^{C) is used instead. These quantities can be calculated by extending the derivation 
of the mean magnetisation for a g-regular graph presented in Ref. [l]. Specifically, the results of that 
paper are that any spin in the system will be in state F with probability, 

P = Fo{P*,H{t)) , (1) 

where H{t) is assumed to vary adiabatically slowly. For increasing H(t), P* is the smallest solution to 
the self-consistent equation, 

P* = Fi{P*,H{t)) . (2) 
Here, the functions Fm{P*,H{t)) (m = 0, 1, . . . , q), 

F^iP*,H{t)) = Y,r (P*)"(l-P*)'^-"->n, (3) 



n=0 

are defined in terms of the probabilities, 



Pn 



I pmh , (4) 

-J{2n-q)-H 

that the local field at a site is positive when it has n fiipped neighbours. The probability density function 
p{h) is assumed to be positive and characterised by a continuous cumulative function. In this case, pn 
are continuous functions of H obeying < po < . . . < pg < 1, implying that < P* < 1. The special 
case of a rectangular distribution, when p{h) can be zero, has been analysed in Ref. [3], where it was 
shown that the discontinuous behaviour in magnetisation associated with an infinite avalanche occurs 
for any regular graph, including a linear chain. 

Extending the analysis of Ref. [l| we can find the rate of field-induced changes in the following form. 
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Similarly, the probability, P^^ , that a spin becomes unstable when its neighbour flips is given by. 



P^^ =q-\i-p*y 



1 fd FoiP*,H{t)) 
dP* 



Using arguments similar to those presented in Ref. [l], it can be shown that 

P^(C,) = {py-^^ (1 - f dFo{p*,H{t)) Y' 



and, 



P^'iCi) = ^. ^) (P*)'^-^-^^ (1 - P*f 



dH 

dFo{P\H {t))Y' 
dP* 



(6) 



(7) 



(8) 



For adiabatically slowly changing H{t), no neighbour of an unstable site i makes the S^U transition 
while site i is in state U. This means that when spin i flips, each of the d stable neighbours of site 
i will become unstable independently with probability P^^ , resulting in the following distributions of 
numbers of offspring, 



P'i^') = E (S) (P'^'f {l-P'''f~''p'iQ and 



q 



The formulae given by Eqs. ©-([lO]) hnk the zt-RFIM to a BP. 



Ci—^ 



(9) 
(10) 



II. DERIVATION OF dE[^]/dt 

In this section, we derive the expression for dE[^]/dt used in the main text. It follows from Eqs. ([2]) 
and ([TO]) that. 



dFi{P*,H{t)) 



l-P* 



dFi{P\H{t)) 



dP* 1 - Fi{P*,H{t)) V dP* 

Et=l (g - 1 - n) (P*)" (1 - Pr-'-"- (Pn+l - Pn) 



The derivative of E[^] can be evaluated as, 

dE[^] _ fdE[^]\ fdP*\ ^ pE[g] ^^ |^dg^^ 



dt \ dP* J \ 'it J \ dH J \ dt J 
where the derivative of P* can be found using Eq. ([2]) as. 



dP* 



1 



dFi{P*,H{t))Y^ fdFi{P\H{t))\ fdH\ 



(11) 



(12) 



(13) 



dt V* dP* J \ dH J \dt J ' 

which is discontinuous only at dFi{P* , H{t))/dP* 1". Combining Eqs. ([TT]) . (fT2]) and ([131) "^e obtain 



dE[^ 
dt 



A{t) 



q-l q-1 

5^5^P^(n,,C.)C, 

.0=0 n,=0 



g-1 g-1 



g-1 g-1 



Hi 



in,; 



, O=0ni=0 



Xi=0 n,=0 



where the function, 
Bit) = 



d'F,{P*,H{t))\^ , fdFr{P*,H{t))\ fdF,{P*,H{t)) 



dHdP* 



dH 



dP* 



dH 
~dt 



Bit) , 
(14) 

(15) 



3 



remains finite at all times. In contrast, the function, 

Ait) = (1 - E[e])- (p* (1 - (^^^^^^^) (^) , (16) 

diverges, A{t) — +oo, as E[^] — )• 1~ because all four factors in Eq. dTUD are positive (see Eq. The 
probability distribution P^{ni,(i) in Eq. ()14p is defined by, 

P^imXi) = P^(C.)^''''5c„,-i-n. + P^(n.)(l - P'^)'5c.,o , (17) 

where, 

PV.) = ~ ^) (^*)"" (1 - Py~''^ (1 - Pn.+l) [(1 - P*) (1 - P'")] ■ (18) 

Here, the probability distribution P^(nj) describes the probability that spin i, which remains stable 
when one of its neighbours change from U to F (i.e. reproduction does not occur at site i, indicated by 
R), has rii flipped neighbours in addition to the one which changed state. The product of probabilities 
pR(Ci)P^u in Eq. ([IT]) is the probability that a spin i becomes unstable and has Q stable neighbours when 
one of its neighbours moves from state U^F. The Kronecker-(5's in Eq. (|17p ensure that Q = q — 1 — n 
when spin i becomes unstable and when spin i remains stable. The probability distribution P^{ni, Q) 
is therefore the joint distribution of rii and Q for any stable spin for which a neighbour changes state 
from U^F, meaning that the term in the square brackets in Eq. (jl4p is an expression for the covariance, 

Cov(C., ni) = K^iQrii] - E^[Ci]K^[ni] . (19) 

used in the main text. Comparison of Eqs. (jlOp and ()17p reveals the formula E[^] = E^[^j], also used in 
the main text. 



III. GENERATING FUNCTION APPROACH TO DYNAMICS 

In this section, we derive equations for the spectrum, pulse shape and distribution of durations of 
avalanches, using the generating function formalism. 

During a single avalanche, corresponding to a particular life-period of the BP process with immi- 
gration, several spins can flip. Let us consider one of these avalanches consisting of flipped spins {sj}. 
Within the continuous-time dynamics, these spins Si will flip at different times ti > 0, with the flrst flip 
occurring at time t = and the flnal one at time t = T. The quantity of experimental interest is the 
voltage on a coil around the sample 0], V{t), which is proportional to the derivative of the magnetisation 
and can be deflned as V(t) = Vo6{t — ti) with a constant of proportionality Vq. 

Let us assume that a single spin flips spontaneously at t = to initiate the avalanche (only one 
fleld-induced change). When this spin flips, a random number, of its neighbours become unstable, 
which is distributed according to P^{^^). Assume that spin i in the avalanche flips at time tj. This 
spin will also have a random number, ^, of neighbours in the down-state which become unstable and 
this number is distributed according to P^{^). Because H{t) and thus P* vary adiabatically slowly 
with time, the values of P^{£,^) and P^{£,) remain unchanged during an avalanche. Accordingly, the 
spin-flip dynamics for any avalanche in the zt-RFIM can be mapped to a simpler form of BP, known 
as a continuous-time Galton- Watson (GW) process 0]. Therefore, we apply the generating function 
formalism developed for the GW process in order to describe the dynamical behaviour of the zt-RFIM. 
Note, however, that in general both P^{C^) and P^(^) vary between different avalanches. In particular, 
this implies that the first moment of P^(^) given by E[^] = E^[Cj] changes during the course of complete 
magnetisation reversal (cf. Fig. 2 of the main text). 
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The generating functions corresponding to P^{^^) and -P^(^) can be defined as, 

q 

a{x) = K[x^'] = ^ P\C^)x^' , (20) 



n=0 



and, 



/(x)=E[x«] = 5]pi^(e)x« , (21) 

where P\^^) and P^(C) are given by Eqs. ^ and (jlOp . respectively. The number N'^ of unstable spins 
in the avalanche at time t is randomly distributed according to P{N^ ,t), and a generating function can 
also be written for this quantity, 

oo 

Gt(x) =E[x^"] = ^ P(7VU,t)x^'' . (22) 

NV=o 

Each spin i unstable at a time t will subsequently flip at time ti > t. This flip can cause its neighbours 
to become unstable, which will subsequently flip causing further spins to become unstable etc. In this 
way, there will be a total of Xi^t' unstable spins at a time t + 1' emerging from a single initial unstable 
spin i at time t. The value of Xi^t' is randomly distributed with probability P{Xi^t'), identical for all i 
and independent of the time t. The corresponding generating function is 

oo 

/t,(x) =E[x^M'] = ^ • (23) 

Each of the unstable spins i at time t = leads to Xi^t unstable spins at a subsequent time t > 
where each of the Xi^t is distributed in the same way as Xi^f. Thus the total number of unstable spins 
at a subsequent time t > is (t) = ^f^i Xi^t and its generating function can be written as, 

Gt{x) = E[x^"W] = E[x^-i^>>*] = E[{E[x^^'']f] = a{ft{x)) . (24) 

Using techniques previously applied to the GW process |4| the generating function ft (x) can be found 
to obey, 

^Al^ = r[/(/,(x))-/i(x)], /o(x) = x, (25) 
which means that it is the solution of the following equation. 



fix') 



dx' = Ft . (26) 



if f{x) — X ^ 0, and ft{x) = x if 



fix) - X = (27) 



The function /(x) takes the value 1 for x = 1 and is convex in the interval [0, 1] since it is a polynomial 
with positive coefficients, P^(^) (see Eq. 1^2)) ■ Therefore, Eq. (glD can have either one non-degenerate 
solution, X = 1 (stable, dxfi^) < 1), a single degenerate solution at x = 1 (stable for perturbations 
described by x < 1, 9^/(1) = 1, d^xfi^) > 0) or two solutions, x = Q < 1 (stable, dxfiQ) < 1) 
and X = 1 (unstable, dxfil) > 1). This implies that there could be three possible regimes for the 
non-stationary solutions of /t(x) given by Eq. (j26|) . sub-critical, critical and super-critical, respectively, 
which correspond to popping, crackling and snapping behaviour in the zt-RFIM (see Fig. 31]) • It should 
be mentioned that the critical regime for the BP occurs not only at the critical point, A = Ac and 
H = He, of the zt-RFIM. In fact, the BP enters the critical regime as the external field approaches the 
coercive field at which snapping occurs, //cocr(A)5 for any disorder A < Ac. The function /(x) — x is 




Fig. SI. Denominator of integral in Eq. ()26p vs x for high disorder, A = 3.0 and H = 1.0 (red-dashed 
hne, sub-critical regime of BP), at the critical point of the zt-RFIM characterised by Ac = 1.78215 
and He = 1.0 (black solid line, critical regime of BP), and immediately before snapping at low disorder 
A = 1.0 and the coercive field i?coer(A = 1.0) = 1.44754 (blue dot-dashed line, critical regime of BP). 
The green double-dot dashed line would correspond to a hypothetical super-critical regime of the BP 
process but cannot be achieved for any values of H or A. 



given by Eq. (j2ip and it is non-negative for < x < 1 for all values of if or A (see Fig. ^T]) . This means 
that the BP cannot be in the super-critical regime. Indeed, any small increase in H from H = //cocr(A) 
will cause the BP to enter the super-critical regime but an infinite avalanche (corresponding to snapping) 
will occur immediately, thus preventing the system from remaining in that regime. 

Using Eq. (j26|) . we analyse the behaviour of the function ft{x) in the sub-critical and critical regimes 
of the BP. It follows from the fact that the integrand in Eq. (|26p is positive and diverges at x' = 1 
that the value of the generating function ft{x) lies in the interval x < ft{x) < 1 for any value of x in 
the interval < x < 1 and t > 0. For such values of x, the value of ft{x) can be found numerically 
in the following way. First, the solution, P*, to the self-consistent equation ([2|) should be found. The 
coefficients of the polynomial function f{x) are then given by Eq. (|lUp . The integral in Eq. (j26p can be 
evaluated in terms of the roots of the polynomial f{x) — x and their residues. Values of ft{x) can then 
be found by inverting Eq. ()26p numerically, while derivatives of ft{x) with respect to both x and t can 
be written in closed forms in terms of ft{x). 

When ft{x) — 1, corresponding either to x = 1 or large values of t, it is possible to write /t(x) and 
its derivatives in terms of P* by expanding the integrand of Eq. ()26p in a Taylor series around the point 
of divergence at x' = 1, i.e. 



'7rT\ — 7^^' - 

'=x f{x') - x' 



1 d^^f(l) 



dx' 



(9,/(l) - 1) (x' - 1) 2(a,/(l)-l)^ 
E[C - 1]-^ (ln(l - ft{x)) - ln(l - x)) - 2-^¥.m " " ^~\ft{x) - x) , 

(28) 

where Eq. (fTOj) can be used to calculate the coefficients of the polynomial /(x), giving, 

d.m - 1 = E[e - 1] = ^ - 1 (29) 

d.J{l) = Em - 1)] = (1 - Q-') ■ (30) 
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For X — )• 1, Eq. ()28p gives for the sub-critical regime of the BP, 

ft{l) = 1 
d,ft{l) = exp [TtE[C-l]] 

d..ftil) = ^^^l^I^^^ [1 - exp [nE[C - 1]]] exp [nE[C - 1]] , (31) 
while for large t ^T~^, 

dtft{x) oc exp [rtE[^ - 1]] . (32) 

When the BP is in the critical regime, — 1] = and the expansion at large t ^ T^^ reads instead 
as. 

Using Eqs. ([20]) - ([33]) . the distribution of avalanche durations, the avalanche pulse shape and the power 
spectrum of avalanches can be found. 

The distribution of the duration of avalanches, pT, can be found using Gt{0) = P{N'^ = 0), where 
P{N^ = 0) is the probability that an avalanche is extinct at time T. The distribution of avalanche 
durations is the derivative of this quantity, 

PT = = dMhrndTMo) . (34) 

As explained above, the value of /t(0) and dTfriO) can be found by solving Eq. (j26|) numerically. The 
coefficients of the polynomial a(x) are given by Eq. (j20|) . When the zt-RFIM is away from its critical 
point and is not close to snapping, the BP describing the avalanches is in its sub-critical regime. Thus, 
in the limit of T ^ F^^, px decays exponentially according to, 

PT oc exp[-FTE[l - ^]] . (35) 

Conversely, at criticality in the zt-RFIM, the BP is also in its critical regime, meaning that at large 
times, pt obeys a power law, i.e. 

E[C(e-l)] ^ ' 

for T ^ r~^. Eq. ()34p along with the limiting cases given by Eqs. ()35p and ()36p represent the results 
discussed in the main text. 

The pulse shape can be defined as, 

{V{t))T = (vo Y.S{t-ti)^ , for < t < T (37) 

where the mean is taken over all avalanches with a duration in a small range of values [T, T + 5T]. This 
definition is consistent with those used in discrete time formulation 0] and experimental setups 0]. 
In order to average the pulse shape over all avalanches of fixed duration T, we need to calculate the 
probability that a spin flips within a time interval [t,t + 6t], given that an avalanche goes extinct at 
time T. Let N^{t) be the number of unstable spins at time t. The probability that one of these spins 
flips in an interval [t, t + 6t] (where 5t — > 0) is given by, 

Prob(Flip in [t,t + 6t]\N^ {t)) = TN^ {t)6t + 6{t)6t , (38) 

where the (5-function accounts for the first spin flip which occurs at time t = by definition. If a spin 
does flip it will destabilise ^ new spins and then there will be a total of N^{t) + ^ — 1 unstable spins in 
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the system. After this event, each unstable spin causes a sub-avalanche of spin flips which goes extinct 
before time T with probability /T-t(O). The probability that all of these sub- avalanches will be extinct 
before time T is therefore, 

Prob(Ext < T|Flip in [t, t + 5t] n (t) n = [/T-t(0)]^''^*^+^"^ . (39) 

The probability that the last of these avalanches goes extinct in an interval [T, T + 6T] (where 5T 0) 
is, 

Prob(Ext in [T, T + JT] [Flip in [t,t + 6t] D (t) n = ST— [/T-t(0)]^"^*^+^"^ . (40) 

Combining Eqs. ()38p and ()40p . the probability that a spin flips in the interval [t,t + 5t] given that the 
avalanche goes extinct in an interval [T,T + 6T] can be calculated as, 

Prob(Flip in [t,t + 6t]\Ext in [T,T + 6T]) = pr'^ (vN'^it)-^ [/T-t(0)]'^''(*)+«-^^ 5t + 6{t)6t , (41) 

where the angular brackets indicate an average over the distributions P{N'^,t) and P(0 of N'^it) and 
^, respectively (see Eqs. (f2T]) and ([22]) ). The mean pulse shape is Vq multiplied by the rate at which 
spins flip, i.e. {V{t))T is given by, 



{V{t))T = yoProb(Flip in [t, t + 6t]\Ext in [T, T + 5T]){6t) 



-1 

where the delta-function 6{T — t) accounts for the final spin flip which occurs at t = T, making the 
above equation valid for < i < T. In order to perform the average over the distributions of N'^ and 
Eq. ()42p can be written in terms of the generating functions f{x) and o(x) defined above, 

{Vit))T = yoTpT-'-^ fE[(x/T-i(0))^''W]E[/r_t(0)«] {fT-tiO))-'] , + Vo6{t) + Vo6{T - t) 



x=l 



i^rpT-'^l^ [GtixfT^tmfih^tm ifT-mr^] ^ + vo5{t) + vod{T - 1) 



x=l 



+Vord^f{fT-tm + VoS{t) + Vo6{T - t) , (43) 

where we have simplified the second line in the above expression using Eqs. ()24p and p4p and the 
identities /t(/T-i(0)) = /r(0) and drfTiO) = d^ft{fT-t{0))dTfT-t{0). In order to evaluate {V{t))T 
using Eq. ([ISjl . the coefficients of the polynomials f{x) and a{x) should be found employing Eqs. ([20]) 
and (f2T]l . Then the values of fx-ti^)-, f{fT-t{0)) and their derivatives along with the derivatives of 
a(/T(0)) and ft{fT-t{0)) can be found by solving Eq. as outlined above (below Eq. ([26])). Eq. (HHD 
gives the final expression for the pulse shape given in the main text. The graphs for {V{t))T at the critical 
point presented in the main text are asymmetric for sufficiently small durations T < 10^, preventing 
their scaling collapse as can be seen from Fig. [2l However, for large T the plots collapse onto a single 
scaling function, thus confirming the scaling hypothesis for this system suggested in Ref. 0]. 
The spectrum of the avalanche as defined in the main text can be written as, 

S{u) = { I e''''V{t)dt ^ 

e^t-t') (V{t)V{t'))dtdt' , (44) 

where the angular brackets refer to averaging over all avalanches. The first spin flips deterministically 
at t = and must be treated separately from the other stochastically flipping spins. Therefore, we 
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Fig. S2. Scaling function f{t/T) in the scaling hypothesis for the pulse shape, {V{t))T = T°'f{t/T) 
with a = 1, suggested in Ref . [a, |7| • The pulse shape was calculated for a 4-regular graph with normal 
disorder at criticality (A = Ac, H = 1.0) using Eq. (|43p . Different curves represent different avalanche 
durations as marked in the figure legend, with T = Vq = 1. 



define U{t) = Vq X^j-^i 6{t — ti) as the voltage created by spin flips except the first one, then {V {t)V {t')) 
contains four terms as follows, 

{V{t)V{t')) = Vi5{mt') + V,5{t') {U(t)) + Vo5{t) {U{t')) + {U{t)U{t')) . (45) 

The mean voltage {U{t)) is equal to Vq multiplied by the probability that a spin flips in a time interval 
[t,t + 5t], 

{U{t)) = FoProb(Flip in [t, t + 6t]){Sty^ . (46) 
The mean {U{t)U{t')) contains two terms corresponding to the cases t = t' and t ^ t' , i.e. 

{U{t)U{t')) = yo^Prob(Flip in [t,t + 5t])5{t - t'){6ty^ 

+ yo^Prob(Flip in [t, t + 6t]r\ Flip in [t', t' + 6t\){5ty^ . (47) 

In Eqs. dMI) and ([37|), Prob(Fhp in [t,t + 5t\) is given by, 

Prob(Flip in [t, t + 6t]) = E[r6tN^{t)] = rd^Gt{l)5t , (48) 

and the probability, Prob(Flip in [t,t + 5t] D Flip in [t' , t' + 5t]), that a spin flips in the interval [t,t + 6t] 
and another spin flips in the interval [t',t' + 6t], can be found as, 

Prob(Fnp in [t, t + 6t] n Flip in [t' + 5t]) 
= T\5tf (E[A^U(Wn)(A^U(tmin) - 1)] + E[iVU(t^in)]E[C]) nX\t-t'\] 

= T\5tf[d.,,Gt^,^{i) + a,.Gt„,„(i)a,/(i)]a,/|,_,,|(i) , (49) 

where tmin = min(t,t'). The factor T5tN^ {train), in the first line of Eq. (09]), refers to the probability 
that a spin flips at the earlier of the two times t and t' . After this spin flips, with ^ offspring, there 
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are A^^(tmin) + ^ — 1 unstable spins. Each of these unstable spins gives rise to j unstable spins 

at the later of the two times t and t' , meaning that probability of a spin flip at the later time is, 

Substituting Eqs. (jl5]) -(j39 ]) back into Eq. results in Eq. (3) given in the main text, i.e. 

(f'OO f'OO 
l + 2Td:rail) exp[-nE[l -^]]cos(a;t)dt + ra^a(l) / exp [-nE[l - ^]] 
Jt=o Jt=o 

/ {d,,Gt{l) + d,Gt{l)d^f{l)) / exp [-rnE[l - C]] COS (uju) dudt 

Jt=0 Ju=0 / 



S^ 



E[C - 1]2 +a;2r 



2r-2 



where, 



and, 



So = Vo' 



1 + 



Si = vi 



E[i - e] 



(50) 



(51) 



(52) 



Here, the terms E[^^] and E[((^^)^] are the moments of distributed according to P^{^^). These quantities 
are given by. 



a,a(l) = E[Cl 



5.xa(l)+a,a(l) =E[(e')'] = {l-q 



dP* \dH 



dH 



dFo 
dP* 



dF2\ f dFi 



dH \dH 



-1 



(53) 



Near the critical point, the term E[l — ^] goes to zero according to a power law, E[l-^] oc - Fc|^/^[&], 
while all other moments of ^ and in Eq. ()50p are constant to leading order. This means that near 
criticality and for small values of a;, the spectrum is given by. 



-2/3 



A\H - Hc\'^/^ +uj^T-^ 



(54) 



so that at criticality, the spectrum can be written S{u}) oc Si{H)(jj 2. Integration of the spectrum over 
the whole range of H removes the divergence for small \H — Hc\, so that the integrated spectrum. 



5'mt(^^) 



S{uj,H{t))R^dt 



(55) 



t=— 00 



obeys a power law, 5int(w) oc a;"'^/^, at critical disorder A = Ac. This result follows by substituting 
Eq. ()54p into the integrand in Eq. (j55p . In fact, experimental measurements of the spectrum [3,0] involve 
a sum across avalanches near the critical point, and could therefore follow the integrated avalanche 
spectrum. 
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